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ABSTRACT 



For a subclass of polar-cap models based on electromagnetic cascades induced by 
curvature radiation (CR) we calculate broad-band high-energy spectra of pulsed emis- 
sion expected for classical and millisecond pulsars. The spectra are a combination of 
curvature and synchrotron components. The spectrum of curvature component breaks 
at 150 MeV, and neither its slope nor level below this energy are compatible with 
phase-averaged spectra of pulsed X-ray emission inferred from observations. Spectral 
properties in the combined energy range of ROSAT and ASCA (0.1 - 10 keV) depend 
upon the location of cyclotron turnover energy e ct = K-^-/ sintp in the synchrotron 
component. Unlike in outer-gap models, the available range of pitch angles ip is rather 
narrow and confined to low values. For classical pulsars, a gradual turnover begins al- 
ready at ~ 1 MeV, and the level of the synchrotron spectrum decreases. At ~ 10 keV 
the curvature component eventually takes over, but with photon index a = 2/3, in 
disagreement with observations. 

For millisecond pulsars, the X-ray spectra are dominated by synchrotron component 
with a ~ 1.5, and a sharp turnover into a ~ —1 at e c t ~ 100 eV. 
Relations of pulsed luminosity Lx to spin-down luminosity L s d are presented for clas- 
sical and millisecond pulsars. We conclude that spectral properties and fluxes of pulsed 
non-thermal X-ray emission of some objects, like the Crab or the millisecond pulsar 
B1821-24, pose a challenge to the subclass of polar-cap models based on curvature 
and synchrotron radiation alone. 
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1 INTRODUCTION 

According to recent results from R OS AT and ASCA (Becker 
& Triimper 1997 and Saito 1998, respectively) luminosities 
of pulsed and unpulsed components in the X-ray emission 
from pulsars are related to spin-down luminosities L s[ j, sug- 
gesting thus a rotation-powered origin of X-rays. The un- 
pulsed emission is usually interpreted as synchrotron radi- 
ation from an unresolved nebula surrounding the pulsar. 
A magnetospheric wind of ultrarelativistic particles (with 
Lorentz factors about 10 7 ) will lead to X-ray emission pro- 
vided the nebula contains magnetic field of the order of 
~ 10~ 4 G (e.g. Manning & Willmore (1994) estimate the 
magnetic field strength of 0.17 x 1(U 4 G for the hypotheti- 
cal nebula around B0950+08). 

Pulsed components are of particular interest, since they 
are the direct signatures of non-thermal processes within a 
magnetosphere. Very likely, pulsed non-thermal X-ray emis- 
sion is nothing but a low-energy tail of gamma-ray emission, 



which in most models is a superposition of curvature (CR) 
and synchrotron (SR) emission. In a subclass of polar cap 
models developed by Dermer & Sturner (1994) gamma-ray 
photons have different origin - they are thermal UV/X-ray 
photons subject to Lorentz-boost by inverse Compton scat- 
tering (IC) with beam particles. 

In polar-cap scenarios pulses in gamma-rays should be 
accompanied by pulses in X-rays, with similar shapes, and in 
phase. With such criteria the Crab pulsar would be a model 
example. However, one has to show first that its contribu- 
tions to X-rays and gamma-rays are in right proportions, i.e. 
in (at least qualitative) accord with model predictions. 

Phase-averaged broad-band {X-y) spectra of six, out of 
seven, gamma-ray pulsars are relatively steep, with photon 
index ax 7 decreasing with pulsar's characteristic age (see 
Thompson et al. 1997 for recent review): 
For the Crab ax-, ~ 2, and therefore, luminosities per loga- 
rithmic energy bandwidth in X-rays and in gamma-rays are 
comparable. For B1055-52, which is the oldest gamma-ray 
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pulsar, ax-y ~ 1.6, and its pulsed X-ray luminosity Lx is 
roughly one per cent of its L 7 . 

[Note: To make this comparison we followed the convention 
used in gamma-rays (fi 7 = 1 sr is the solid angle of gamma- 
ray emission) taking Qx = 1 sr for the solid angle of X-ray 
emission.] 

These observed spectral properties have been quite suc- 
cessfully accounted for in a newly proposed model of 'thick 
outer gap' by Zhang & Cheng (1997). Unlike polar-cap mod- 
els, outer-gap models [see also Wang et al. (1997)] may take 
advantage of the full range of pitch angles, which in turn 
makes it possible to trigger electromagnetic cascades even 
for relatively soft gamma-ray photons 1 MeV) on mag- 
netic field lines (also on closed ones). The energy distribution 
function for e* pairs produced in the cascades may easily be- 
come steep enough (N('j±) oc 7 ± ) over a wide range of en- 
ergy, for subsequent synchrotron spectrum (ctx-y = 2) to ac- 
count for the observed soft gamma-rays and X-rays. In par- 
ticular, Cheng et al. (1998) managed to reconstruct simple 
empirical relation between Lx and L s( j , which reads (assum- 
ing fix = 4tt sr) L x ^ 10" 3 L sd (Becker & Trimmer 1997). 
Moreover, cascades propagating starward are subject to sub- 
stantial lateral spread, leading thus to rather weak X-ray 
modulations. 

The aim of this paper is to present general features 
of broad-band X7 spectra expected for the subclass of po- 
lar cap models, with X7 emission due to CR and SR (e.g. 
Daugherty & Harding 1982, Daugherty & Harding 1996) and 
to show that luminosity of pulsed X-rays observed for some 
pulsars is too high compared to L 7 (or to _L s( j if there is 
no information about gamma-rays) to be understood within 
polar-cap models unless some unorthodox assumptions are 
accepted. Some of these features are only weakly model- 
dependent and may, therefore, play a decisive role in asess- 
ing validity of polar-cap models. We were motivated by the 
fact, that neither proponents nor oponents of the polar-cap 
models seem to acknowledge a real challenge posed to these 
models by observed properties of non-thermal, pulsed X-ray 
emission. It is this emission which is more promising as a 
potential discriminator between rivalring classes (polar-cap, 
outer-gap) of models than gamma-ray emission itself. 
In Section 2 we present model spectra of CR and SR emis- 
sion per parent (beam) particle, expected in pulsars with 
high (~ 10 12 G) and low (~ 10 9 G) magnetic fields. Section 3 
compares Lx expected within energy band 0.1 keV— 10 keV 
with observations for a wide range of spin-down luminosities 
L s d- Section 4 contains discussion and conclusions. 



2 SPECTRAL PROPERTIES OF CURVATURE 
AND SYNCHROTRON RADIATION 

The energy distribution JV e (7) of particles is governed by 
their injection rate Qe{"i) (also called a source function), 
cooling rate 7, and a characteristic time scale of their escape 
t csc , via steady-state kinetic equation 



iVe( 7 ) 



|-(| 7 |jVe) =JV e /*e 



Qc 



(1) 



where 7 is particle energy in units of m c c 2 . 

The solution of eq.nl yields two asymptotic relations 



1 p 
TtT Jj 

Cjc ' tet 



(2) 



rmax Q c (7)d 7 , |7l/7>t CS c; 

171/7 < ^se- 
lf the solution to eq|l] is a power-law function with index p, 
N c oc 7 _p , the subsequent photon spectrum of SR or CR 
will have a form 

p+i 



N„(e) oc e 



(3) 



(but the value of p must not be too low) where q = 2 for 
SR, or q = 3 for CR. 

[Note: Throughout this section we use a short notation for 
electron and photon distribution functions: iV e (7) = dNe/d-y 
and iV„(e) = dN v /de, respectively.] 



2.1 Curvature Radiation 

First, let us consider spectral component due to CR. If the 
source function Qc{"i) of beam particles is monoenergetic 



= (7) = Co<5( 7 o - 7), 



(4) 



and redistribution over the energy space is governed by the 
CR cooling (| 7 1 = 7cr j oc 7 4 ), then from the first relation of 
eq.^ it follows that N e oc 7 -4 (i.e. p = 4) for 7 < 70. The 
lower limit (let us denote it as 7break) to this distribution 
is determined by the condition, that a cooling time-scale 
due to CR, t C r = 7/| 7c r|, is shorter than t csc . According to 
eq.^, the photon spectrum of CR becomes iV cr (e) oc . 
It remains unchanged down to photon energy ebrcak, which 
may be approximated with a characteristic energy of CR 

3 7 3 

Ccrit = - cTl , (5) 

2 Per 

(where p cr is a local radius of curvature) taken for 7break- 
Below ebrcak, the photon spectrum follows the low-energy 
tail of CR, i.e. AT cr (e) oc e~3. In order to find ebrcak we 
should first estimate the time scale t eac . We assume that 
t C sc ~ Pcr/c. The CR cooling rate for a single particle of 
energy 7 is 



2 fix* 

3 mc pi 

The condition t CI = t CS c becomes then 



7break 
Per 



3 mc 



and from eqj^ it follows that 



_ 9 c_ 

^brcak — . " 

4 r 



150 MeV, 



(6) 



(7) 



(8) 



where ro is the classical electron radius. Note, that the pho- 
ton energy tbreak at which the spectral break occurs does 
not depend on any pulsar parameters, as long as our esti- 
mation of tcsc is accurate. In particular, it does not depend 
on magnetic field structure - should it be pure dipole with 
Per <; 10 8 \/P cm, or dominated by high-order multipolar 
components with p cr ~ 10 6 cm. 

Fig.l shows numerically calculated shapes of the CR 
spectrum produced by beam particles injected at the outer 
rim of a canonical polar cap with a dipolar magnetic field. 
The magnetic strength B pc at the polar cap is 10 9 G (upper 
panel), and 10 12 G (lower panel). The period of rotation P 
is 0.003 s and 0.06 s, respectively. For the initial energy Eo 
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we choose 1.1 x 10 7 MeV and 6.7 x 10 6 MeV, respectively 
(see section 3 for explanation). The dotted line corresponds 
to the unabsorbed spectrum, and the solid line is the spec- 
trum with magnetic absorption. The spectral break (where 
the power-law spectrum changes its slope by one power of 
e) is visible in both panels, and its location at ~ 150 MeV 
agrees quite well with the analytical estimate (eq.^) . A word 
of comment on high-energy cutoffs seems appropriate here, 
though it is not linked directly to the subject of this paper: 
The high-energy cutoff for the case of 10 12 G does not ex- 
ceed 10 GeV, whereas the cutoff for the low-B case reaches 
~ 0.1 TeV (the energy range accessible with ground-based 
Cherenkov techniques). The explanation is twofold: (1) for 
low values of B the magnetic absorption is weak; (2) very 
short periods P infer small dipolar curvature radii, and a 
characteristic photon energy (eqjEJ) taken for 70 = Eo/m c c 2 
increases. 

With the spectral break at tbrcak ~ 150 MeV, the cur- 
vature radiation becomes energetically unimportant below 
this energy. Such objects would be dim in X-rays, which ap- 
parently is not the case. Moreover, the slope —2/3 in the 
X-ray energy range is in clear conflict with spectral analysis 
of X-ray data (Becker & Triimper 1997). 

A way out may be offered by including synchrotron ra- 
diation (next subsection) and/or by relaxing the assumption 
about the monoenergetic source function of beam particles 
(eq^J). Suppose, for example, that we want the slope of —5/3 
to extend down to energy e x in the X-ray range: N v (e) oc e~ 3 
for e <; e x . According to eq.^ (lower relation) we now have 
the energy distribution N c (-y) — Q c ■ t asc . This, in turn, re- 
quires that the source function of beam particles is 



Q c oc 7 4 , for 7 i> -y x , 



(9) 



a very steep function of 7 (for whatever the reason), where 
j x is linked to e x via eqM To find j x we solve the equation 
e x = e C rit(7a:) for *y x : For dipolar magnetic fields, the smallest 
curvature radii p CI are those at outer rim of the polar cap. 
They may be approximated with p CI ~ \/r Rl, where r is 
the radial coordinate of the particle, and Rl = cP/2ti is the 
pulsar's light cylinder. Within the region of r ~ 2Rns we 
have therefore t CBC « 10~ 2 ' 5 \/P s. From eq.|| we obtain 

e x ^0.3( 7a; /10 5 ) 3 p- - 5 keV. (10) 

For a pulsar with P — 0.1 s, the source function Qeipi) oc 7~ 4 
would have to extend down to j x ~ 10 5 for e x to reach 
~ 1 keV. This value of y x is about 55 times lower than the 
value of 7brcak (eq.|) corresponding to e brC ak (eqj|) . 



2.2 Synchrotron Radiation 

Since CR is not a promising candidate for hard X-rays, the 
standard explanation involves SR, produced by e* pairs cre- 
ated in the process of magnetic absorption of high-energy 
CR photons. Total contribution of the SR component to 
the overall (SR+CR) energy of radiation per beam parti- 
cle depends on number of created pairs n±. One may ex- 
pect, therefore, that at least for pulsars with strong mag- 
netic fields (~ 10 12 G), the energy contained in SR should 
be comparable or exceed the energy contained in CR. How- 
ever, spectral properties of the SR component depend not 
only on local values of B, but also on the available range of 



pitch angles t[> between magnetic field lines and the direction 
of propagation of created pairs. The most popular assump- 
tion - about an isotropic distribution - may be fully justified 
in the case of outer-gap models (Zhang & Cheng 1997), but 
is not obvious by any means for polar-cap scenarios. 

Below, we'll present numerical calculations of syn- 
chrotron spectra due to e ± -pairs which are assumed (after 
Daugherty & Harding 1982) to be directed along the di- 
rection of propagation of their parent CR-photons at the 
moment of the creation. This assumption, along with nar- 
row opening angles of magnetic field lines available for beam 
particles, results in a confinement of pitch angles to a nar- 
row range of low values. For a pulsar with a dipolar mag- 
netic field and rotating with period P, we expect then 
sinip < 0.01/ \[P as a first approximation. 

The character of the source function Q± of e ± -pairs 
depends primarily on the reachness of the cascades. In the 
case of low magnetic fields (~ 10 9 G) only one generation of 
e ± -pairs is produced, and their formation region is narrow. 
For these reasons Q± is almost monoenergetic (see Fig.l, 
upper panel) , and basic spectral properties of SR calculated 
for a canonical millisecond pulsar are easily understood with 
well known analytical considerations: 

Let us denote by 711 the particle Lorentz factor along the 
local magnetic field line. We will ignore possible curvature of 
the line, and assume that 711 remains constant. The particle 
energy (in units of m c c 2 ) to be emitted in bursts of SR 
corresponds then to Lorentz factor j± of gyration. The total 
particle energy is 7 = "/± 711. As long as 711 3> 1, the pitch 
angle of the particle is determined by siwip ~ 7i| • 

The rate of SR cooling is 



2 r 2 a 2 

7sr = -- B 7 ± 

3 m c c 



2 ^0 D 2 • 2 , 2 
B sm tp 7 . 

3 m c c 



(11) 



Unlike in the case of CR, the SR cooling rate is enormous 
and affects the energy distribution of pairs until 7x ~ 1 
where the synchrotron approximation breaks (O'Dell & Sar- 
tori 1970). Therefore, from eqsj|, [n] and H, we'll get pho- 
ton spectrum of SR with well known single power-law shape 
iV sr (e) oc e~2. Main contribution to the spectrum at energy 
e comes from particles with 7 for which e = e S r(7), where 



- cn 7 simp, 

2 m e c 



(12) 



is a critical photon energy in SR. The spectrum spreads 
between a high-energy cutoff e sr (7max) and a low-energy 
turnover e ct determined by the condition 7^ ~ 1 (O'Dell 
& Sartori 1970): 



_ 3 eB 1 

e rt = e S r(7 = 7|| ) = ^ ch 



m e c sin tp 



(13) 



Below e ct , the spectrum flattens, and may be described 
asymptotically as N BT (e) oc e +1 . It is built up by con- 
tributions from low-energy tails emitted by particles with 
71 > 1, and each low-energy tail is assumed to cut off at 
local gyrofrequency, which in the reference frame comoving 
with the center of gyration is ujb — 



The numerical example presented in the upper panel of 
Fig.l reveals all these 'classical textbook' features so clearly, 
because (for reasons mentioned before) the source function 
of pairs in the case of a dipolar field with B pc — 10 9 G 
is practically monoenergetic. The turnover energy is located 



© 0000 RAS, MNRAS 000, 000-000 



4 B. Rudak & J. Dyks 



- 









ci 






« 


Q 




iafl 


C 





-5 5 
— i 1 1 1 1 1 1 1 1 1 1 




- 

2 

4 

6 



- 

- -2 
4 

-6 



-5 







log & [MeV] 



Figure 1. The radiation energy spectrum per logarithmic energy bandwidth. The spectrum is normalized to the energy of the parent 
particle, Eq. It consists of two components: curvature and synchrotron. Dotted line is the curvature spectrum before correction for the 
magnetic absorption effects. Solid line connecting filled dots, is the curvature spectrum with the magnetic absorption taken into account. 
Solid line without any symbols overlayed is the synchrotron component. Dot-dashed line indicates superposition of the two. 
In addition to the electromagnetic spectra, we show the spectra of e^ pairs: the line connecting filled squares is for the energy parallel to 
local magnetic field lines (labeled with the line connecting open squares is the initial distribution of energy perpendicular to local 
magnetic lines (labeled with £±). 

The upper panel is for B pc = 10 9 G, and E = 1.08 X 10 7 MeV (see eq.^) which corresponds to P = 3.1 X 10~ 3 s (i.e. L sd = 10 35 erg s _1 ). 
The lower panel is for B pc = 10 12 G, and E = 6.68 X 10 6 MeV which corresponds to P = 5.6 X 10 -2 s (i.e. L sd = 10 36 erg s" 1 ) 



around 0.1 keV, i.e. in the soft range of X-rays. More impor- 
tantly, the SR component dominates over the CR component 
in the entire energy range of X-rays. The lower panel of Fig.l 
shows analogous results obtained for B pc = 10 12 G. The cre- 
ated pairs belong now to two generations, and they are rich 
(~ 400 times more numerous than for 10 9 G). Therefore, 
the SR component is now energetically comparable to CR. 
There is no single, well denned turnover energy e ct anymore. 



The spectrum of the source function Q± is spread over two 
decades in energy, and in consequence the SR spectrum re- 
veals a gradual turnover, which starts already at ~ 1 MeV, 
due to high values of 7y as well as strong local B (see eq.|l3|). 

Our approximate treatment of SR spectral shapes at 
low-energy limit requires a word of comment: As electrons 
go to the ground Landau level, the spectrum reveals its har- 
monic structure and the analytical formula that we used, 
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should be treated with caution. Comparison of Monte Carlo 
spectra (calculated by summing-up the rates of the quan- 
tum transitions), with analytical formulae was carried out 
by Harding & Preece (1987). They found that the analyt- 
ical formula for the low-energy tail of SR with a cutoff at 
cyclotron frequency generally overestimates (but not dra- 
matically) the actual level of photon spectrum. Therefore, 
our results for SR in the context of X-rays should be treated 
rather as upper limits, but this does not change our conclu- 
sions. 



3 RELATIONS BETWEEN X-RAYS AND 
GAMMA-RAYS 

We calculated combined spectra of CR and SR (as shown in 
Fig.l), emitted by a beam particle with initial energy Eq, 



for two canonical values of B„ 



10 a G and 10 i2 G, and for 



a range of periods P in order to cover a full possible range of 
spin-down luminosities L s( j. Formula for the energy Eo was 
adopted from Rudak & Dyks (1998) [RD98]. This energy is 
only a few times higher than threshold energy -E m i n required 
for magnetic pair creation in a dipolar magnetic field: 



E = 2.5 x E n 



with E n 



1.2 x 10 



*(-*-) 

Vio 12 g) 



-1/3 



,1/3 



McV, 



(14) 
(15) 



but at the same time it must never exceed any of the fol- 
lowing limits - i?max and E w (the first restriction is due to 
curvature cooling, the second one is due to potential drop 
across the polar cap). The original motivation for this for- 
mula was to reproduce L 7 for the seven gamma-ray pulsars. 
For 10 9 G, E « 10 7 MeV, while L sd covers the range be- 
tween 10 34 erg s _1 and 10 37 erg s _1 . In practice, the values 
of Eo do not change dramatically over the allowed ranges of 
L sd : For 10 12 G, E starts with 3 x 10 s MeV, increases to 
9.8 x 10 6 MeV, and goes down to 4 x 10 6 MeV, for L sd chang- 
ing from 10 31 erg s _1 to 10 39 erg s -1 . 

What became clear already in the previous section is 
that the range of X-rays is never energetically important. 
Although for 10 9 G the spectrum of SR extends well into soft 
X-ray band and dominates there, its fractional contribution 
to the total radiation energy output is low. It varies between 
0.003 (for L sd = 10 34 erg s" 1 ) and 0.22 (for 10 37 erg s" 1 ). 
For the case of 10 12 G, the energy content of SR becomes 
significant for L s d > 10 33 erg s _1 , but the spectrum of SR is 
now confined to gamma-rays, and it turns over at ^ 1 MeV. 
In either case, the bulk of particle energy converted into 
radiation is concentrated within gamma-rays. 

Gamma-ray luminosity L 7 in RD98 was identified with 
the power of outflowing particles 



^particle; 



= r]E±n±N GJ , 



(16) 



where E± (characteristic energy attained by a fraction 
of secondary particles, rjn±, due to possible acceleration) 
was assumed to be of the order of Eo, and Ngj was the 
Goldreich- Julian rate of outflow of beam particles. The pa- 
rameter r\ — 0.004 was used to match the model with 
the joint EGRET and COMPTEL luminosity inferred for 
B1951+32 (see RD98 for detalis). Here we made more ac- 
curate calculations for L 7 , by replacing the particle energy 
Eq with its radiation energy yield Ej for photon energy 



30 32 34 36 38 




30 32 34 36 38 40 
log L sd [erg s" 1 ] 



Figure 2. Evolution of high-energy, non-thermal luminosity 
across the spin-down luminosity space. The long upper curve con- 
necting rilled dots is the track of gamma-ray luminosity L 1 {> 
100 keV), and the long lower curve connecting filled squares is 
the track of X-ray luminosity Lx(0.1 — 10 keV), both calculated 
according to cq,|l7| for a pulsar with B pc = 10 12 G. [The dotted 



(cq.[u^), which in RD98 was identi- 



line is the track of L partic i es 
fled with Lj.] 

The short upper curve with open circles is for gamma-ray luminos- 
ity L 1 {> 100 keV), and the short lower curve with open squares 
is for X-ray luminosity Lx(0.1 — 10 keV), calculated according to 



eq. 







for B„ 



10 9 G. The short-dashed line marks L = L„, 



The long-dashed line, marking the empirical relation of Becker 
&; Triimper (1997) rewritten for !2j = 1 sr, has been added for 
reference. The fluctuations of L x (0.1- 10 keV) (for 10 12 G only) 
between 10 36 erg s — 1 and 10 37 erg s — 1 in Z/ sd , are of numerical 
origin: in this range of L sd the low-energy tail of the synchrotron 
component crosses the curvature component just around 10 keV, 
and any fluctuations due to Monte Carlo treatment of cutting off 
synchrotron spectra influence the level of the total spectrum. 



e > 100 keV. Similarly, we calculate expected X-ray lumi- 
nosity Lx within (0.1 keV — 10 keV). Accordingly, we take 



L 



particles X 



Eo~ 



and Lx = L 



particles ^ 



Ex_ 
Eo 



(17) 



where ^ and are fractional radiation energy yields per 
particle, calculated by integrating the differential (per loga- 
rithmic energy bandwidth) spectra [which are simply 

1 dN v . . 

Wo'd^e (18) 
- see Fig.l for two examples], over the range e > 100 keV 
and 0.1 keV < e < 10 keV, respectively. 

Extremely low values of relative energy output in X- 
ray photons and gamma-ray photons per particle imply 
equally low values of Lx /L~, (between 10 -6 and 10~ 4 ). Fig. 2 
presents L 1 and Lx calculated according to eqs.[l7j, [lcj , and 
hq as a function of L sd . It shows unambiguously, that some 
pulsars from the list of X-ray and gamma-ray source do not 
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fit the picture: The Crab pulsar, for which the observed ratio 
of Lx /L-y is about 0.47, remains certainly a challenge. Also 
among millisecond pulsars, there is a 'Crab-like' example, 
which does not fit the model. B1821-24 observed with ASCA 
(Saito et al. 1997) reveals clear, double-peaked pulsations 
resembling those of the Crab. Apart of the disagreement be- 
tween the observed spectral slope (~ 1.9) and the predicted 
one (~ 1.5), the predicted ratio of Lx/L 7 ~2x 1CT 4 (we 
took Lx for lsr) would make L 7 to exceed L s d. Or, alter- 
natively, using the upper limit for the gamma-ray luminos- 
ity, L 7 < 0.028 L sd (Nel et al. 1996), the "observed" ratio 
Lx/Lj is > 2.5 x 10~ 3 , in disagreement with predictions. 

However, among pulsars there are also strong sources of 
steady X-ray emission, with no trace of pulsations. B1706- 
44, a strong gamma-ray pulsar detected with EGRET, with 
L 1 = 2.5 x 10 34 erg s _1 (Thompson et al. 1996) shows no 
pulsations in X-rays. Its unpulsed, non-thermal X-ray emis- 
sion has the luminosity Lx = 1.3 x 10 33 erg s™ 1 (Becker et 
al. 1993). Apparently, pulsed X-rays of B1706-44, which we 
predict to reach the luminosity of ~ 3 x 10 erg s~ , arc 
dwarfed by nearby nebular steady emission. 



4 CONCLUSIONS 

Recent results from ROSAT and ASCA suggest that the X- 
ray spectra of pulsars may be a superposition of thermal 
and non-thermal components (see e.g. Zavlin & Pavlov 1997 
for discussion). The existence of these two components is in 
general expected in both classes of rival models (polar-cap 
and outer-gap models; e.g. Arons 1981 and Cheng, Ho & Ru- 
derman 1986, respectively) and the thermal component re- 
sults from polar cap heating by inflowing particles. However, 
clear empirical relations between thermal and non-thermal 
components are yet to be determined. The latest review on 
the status of pulsar' X-ray properties by Becker & Triimper 
(1997) favoures pure non-thermal spectral models (with four 
exceptions of additional thermal components due to initial 
cooling of the star itself). 

Within the framework of polar-cap models with mag- 
netospheric activity induced by curvature radiation of beam 
particles we have calculated broad-band photon spectra of 
non-thermal origin, for both classical and millisecond pul- 
sars. This non-thermal emission is a superposition of curva- 
ture and synchrotron radiation. It is expected to be pulsed, 
like its high-energy extension - the gamma-ray emission. Our 
objective was to estimate spectral properties and the level 
of this emission in the energy domain of X-rays. 

Our calculations show, that for beam particles with mo- 
noenergetic source function, the properties of photon spectra 
in the combined energy range for ROSAT and ASCA (0.1 
- 10 keV) depend primarily on the magnetic field strength 
at the polar cap. Millisecond pulsars, with dipolar magnetic 
fields of the order of 10 9 G, have X-ray spectra dominated 
by synchrotron emission. On the other hand, pulsars with 
B of the order of 10 12 G have X-ray spectra almost exclu- 
sively due to low-energy tail of curvature emission. However, 
the spectrum of the curvature component breaks already at 
150 MeV regardless the pulsar parameters, and neither its 
slope nor level in the X-ray domain is compatible with phase- 
averaged spectra of pulsed X-ray emission inferred from ob- 
servations. 



Detailed properties of the broad-band spectrum, de- 
pend upon the location(s) of cyclotron turnover energy 
£ct = h-^:/ s'mip of the synchrotron emission. Unlike in 
outer-gap models, the available range of pitch angles ip is 
rather narrow and confined to low values (sin?/; < 0.01/\/P). 
We find that for classical pulsars (10 12 G) the upper limit 
for turnover energy e c t occurs already at ~ 1 MeV. Be- 
low this energy, the level of synchrotron spectrum gradually 
decreases, and around ~ 10 keV the curvature component 
takes over. Unless the spectrum of primary beam particles is 
very steep (Nb oc % 4 ) and extends down to Lorentz factor 
7b ~ 10 5 , the X-ray spectrum below ~ 10 keV is a power- 
law with photon index a p h = 2/3. Therefore, the model 
is not able to explain the non-thermal X-ray spectra with 
a p h — 1.5-^2 inferred from observations (Becker & Triimper 
1997). 

The expected luminosity Lx of pulsed component is a neg- 
ligible fraction of gamma-ray luminosity L 7 , and for L sd > 

10 33 erg s _1 an approximate relation 

log L x « -1.5 + 0.83 log L sd (19) 

follows from Fig. 2. The situation is qualitatively different for 
millisecond pulsars. The X-ray spectra are dominated now 
by synchrotron components, with a photon index a p h ~ 3/2, 
which extends down to e c t ~ 100 eV before breaking sharply 
into a p h ~ — 1. The expected Lx is still low with respect to 
L 7 , though not as low as for classical pulsars. For L s d > 

10 34 erg s _1 we find 

log L x « -40.3 + 1 .90 log L sd . (20) 

We conclude that spectral properties and fluxes of 
pulsed non-thermal X-ray emission of some objects, like the 
Crab or the millisecond pulsar B1821-24, pose a real chal- 
lenge to the subclass of polar-cap models based on curvature 
and synchrotron radiation alone. 



ACKNOWLEDGEMENTS 

This work has been financed by the KBN grant 2P03D- 
00911. BR acknowledges discussions with K.S. Cheng, J. 
Gil, and W. Kluzniak. We thank the anonymous referee for 
useful suggestions which helped to clarify the paper. 



REFERENCES 

Arons J., 1981, ApJ, 248, 1099 
Becker W., 1995, Ph.D. Thesis (LMU Munich) 
Becker W., Brazier K.T.S., Triimper J., 1993, A&A, 273, 421 
Becker W., Triimper J., 1997, A&A, 326, 682 
Cheng K.S., Gil J., Zhang L., 1998, ApJ, 493, L35 
Cheng K.S., Ho C, Ruderman M.A., 1986, ApJ, 300, 500 
Daugherty J.K., Harding A.K., 1982, ApJ, 252, 337 
Daugherty J.K., Harding A.K., 1996, ApJ, 458, 278 
Dermer CD., Sturncr S.J., 1994, ApJ, 420, L75 
Harding A.K., Preece, R.T., 1987, ApJ, 319, 939 
Manning R.A., Willmore A.P., 1994, MNRAS, 266, 635 
Nel H.I. et al., 1996, ApJ, 465, 898 
O'Dell S.L., Sartori L., 1970, ApJ, 161, L63 
Rudak B., Dyks J., 1998, MNRAS, 295, 337 [RD98] 
Saito Y., 1998 in Shibata S., Sato M., cds. Proc. Neutron Stars 
and Pulsars. Universal Academy Press, Tokyo, in press 



© 0000 RAS, MNRAS 000, 000-000 



High- energy emission 7 



Saito Y., Kawai N., Kamae T., Shibata S., Dotani T., Kulkarni 

S.R., 1997, ApJ, 477, L37 
Thompson D.J. et al. 1996, ApJ, 465, 385 

Thompson D.J., Harding A.K.H., Hermsen W., Ulmer M.P., 1997, 
in Dermer CD., Strickman M.S., Kurfess, J.D., eds, AIP 
Conf. Proc. 410, 4th Compton Symposium. AIP, New York, 
in press 

Wang F.Y.-H., Halpern J. P., 1997, ApJ, 482, L159 

Wang F.Y.- H.. Ruderman M Halpern J. P., Zhu T., 1997, 

preprint (astro- ph/9711283| ) 
Zhang L., Cheng K.S., 1997, ApJ, 487, 370 
Zhavlin V.E., Pavlov G.G., 1998, A&A, 329, 583 



© 0000 RAS, MNRAS 000, 000-000 



